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In this work we extend recent study of the properties of the dense packing of "superdisks," by 
Y. Jiao, F. H. Stillinger and S. Torquato, Phys. Rev. Lett. 100, 245504 (2008), to the jammed 
state formed by these objects in random sequential adsorption. The superdisks are two-dimensional 
shapes bound by the curves of the form \x\ 2p + \y\ 2p = 1, withp > 0. We use Monte Carlo simulations 
and theoretical arguments to establish that p = 1/2 is a special point at which the jamming density, 
pj(p), has a discontinuous derivative as a function of p. The existence of this point can be also 
argued for by geometrical arguments. 



There has been recent interest @, 3, Q, [j| in 
the problem of geometrical packing and surface deposi- 
tion of noncircular objects in two dimensions (2D). This 
problem is intriguing from the theoretical point of view. 
In addition, it finds applications in studies of design and 
control of prcpatterned surfaces with special properties. 
New capabilities to pattern surfaces at the nanoscale, and 
use particles of nanosizes, have promise for development 
of novel biosensors and detectors, applications in elec- 
tronics @, 0], catalysis @, and optics [§], fioj . 

Recently, an interesting study was reported [l[ of the 
densest possible packing of (oriented) "superdisks" de- 
fined by \x\ 2p + \y\ 2p < 1. These shapes are illus- 
trated in Fig. [T] In particular, numerical evidence for 
< p < 1 (where the p = shapes are defined as a limit 
which yields crosses) suggests [l[ that the point p = 1/2 
separates different closed-packed structures. Note that 
p = 1/2 also separates the convex and concave shapes, 
as shown in Fig. [TJ 

Particle deposition at surfaces is typically irreversible, 
and for a theoretical description of their adsorption one 
can use the random sequential adsorption (RSA) model. 
The RSA model, as well as its various modifications, finds 
app lications and has been extensively studied [ill . 





However, most RSA studies have been carried out for 
spherical (circular) and other simply shaped objects. The 
approach to the jamming density, pj, in RSA processes 
is described by the standard Pomeau [33| and Swendsen 
[HI conjuncture which gives the asymptotic results for 
oriented squares and for disks, which are in a gree ment 
with Monte Carlo (MC) simulation results [35| , |36| (ori- 
ented squares), and [IJ, [3lJ (disks). However, for non- 
oriented squares evidence has been reported [38[ that this 
conjuncture might not work. Moreover, the asymptotic 
behavior of the deposit density for objects with concave 
shapes on continuum substrates has not been studied. 
Studies of RSA of objects with zero area, such as rods, 
circular arcs, etc., have reported interesting features of 
the jamming coverage [39[, [ioj . 

In this work we consider RSA of oriented superdisks in 
two dimensions. We use a grid-type MC algorithm [35[ 
which is particularly suitable for evaluating the density of 
the jammed state, because it efficiently treats deposition 




p = 0.3 



p = 0.5 



0.7 






p=l 



p = 1.3 



p — > oo 



FIG. 1: Superdisks shapes for different values of the deforma- 
tion parameter p. 



in small remaining vacancy areas close to jamming; see 
Fig. [2] Similar to the dense-packing results [l[ , we find 
that p = 1/2 is also a special point for the jammed state 
of RSA. In addition to numerical evidence, this conclu- 
sion will also be substantiated by geometrical arguments. 

A superdisk is a 2D case of the surface of a d- 
dimensional superball pj. A superball is defined as the 
volume of the Euclidean space bounded by the surface 
\x\\ 2p + |x 2 | 2p + ... + \x d \ 2p = 1, where the Carte- 

sian coordinates and p is the deformation parameter. Su- 
peballs have full rotational symmetry only when p = 1 
(when they became hyperspheres) . 

For < p < 1/2 superdisks are concave and for 1/2 < 
p < oo they are convex. The p — > superdisk is reduced 
to cross, the p = 1/2 and p = oo shapes are squares, and 
the p = 1 shape is a circle. 

The reason that we focus on the point p = 1/2 is that, 
with the advent of nanotechnology, and with prolifera- 
tion of ex per iments on deposition of proteins [HI], [Til ]. 
HH > [ill j 12 1 we expect that situations will be realized 



when the particle shapes on the surface, change between 
concave and convex depending on the physical and chem- 
ical conditions of the environment. This might affect 
the asymptotic approach to the jamming coverage (an 
issue that requires a separate detailed study). As demon- 
strated here, the change in the concavity also results in 
a nonanalytic behavior of the jamming coverage, Oj(p), 



2 



at p = 1/2. 

In RSA, a superdisk can be deposited at a surface if 
it does not overlap previously deposited particles. Such 
adsorption is a noncquilibrium process, and therefore the 
deposited particle density does not reach the maximal 
dense packing. Instead, it approaches the jamming den- 
sity, pj(p), at large times. This quantity, the density 
of the deposited particles per unit area, is related to 
the jamming coverage 9j{p) — the fraction of the cov- 
ered area — via 9j(p) = A(p)pj(p), where A(p) is the 
superdisk area. The latter quantity is given by 




where T(x) is the standard gamma function. Since this 
function is analytic near p = 1/2, the behavior of 9j{p) 
and pj(p) at p = 1/2 is easily related. We focus on pj(p), 
because it simplifies some notation below. 

In our MC simulations we used an algorithm originally 
introduced in [35| , which allows to simulate the formation 
of the jamming state, and to estimate pj, using minimal 
computer resources. We used a square system of size 
500D x 500D with periodic boundary conditions, where 
D is the "diameter" of the superdisks along the x and 
y axes, equal 2 in our dimensionless units. Each value 
of p(p) was obtained by averaging over 1000 independent 
runs. The maximum fractional uncertainty in our simu- 
lation was estimated as Apj(p)/ pj(p) — 0.00223. Specif- 
ically, for the squares (p = 1/2) and disks (p = 1) we ob- 
tained the estimates 6j(p = 1/2) = 0.5620 ± 0.0001 and 
9j(p = 1) = 0.5468 ± O.OOOf^which are consistent with 
the values reported in (35[, [3(|, [2l[, (37[. The behavior 
of the jamming density as a function of the deformation 
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FIG. 2: (a) Superdisks (dark shapes) with their exclusion ar- 
eas (lighter shapes). Upper panel: a concave superdisk for 
p = 0.3. Lower panel: a convex superdisk for p = 0.7. The 
dashed lines mark the p — 0.5 squares and their exclusion 
areas, (b) A typical configuration of concave superdisks near 
the jammed state, with at most a single additional superdisk 
deposition possible with its center landing in the central un- 
shaded area. 



parameter is shown in Fig. [3l 

Our data clearly indicate existence of a special point at 
p = 1/2. The p-derivative of the jamming density p,;(p) 
at this point has a discontinuity, similar to that men- 
tioned in [lj for the dense-packing density. In order to 
understand the origin of this behavior in RSA, let us con- 
sider the exclusion area of the superdisks, S(p), defined 
as the area within which it is impossible to deposit an- 
other superdisk's center without overlap, Fig. [2j Unlike 
A(p), the area S(p) markedly changes its p-dependence 
for p above and below the square-shape value of 1/2. It 
is a continuous function of p, but has a discontinuous 
derivative at p = 1/2 (we give the expressions shortly). 
Therefore, on dimensional grounds it is tempting to con- 
jecture that the following relation provides a good quali- 
tative approximation for the superdisk jamming density 
ratio, 

pM „ 3(1/2) m 
P , 7 (l/2) " S{p) 1 ' 

at least near p = 1/2. 
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FIG. 3: Lower panel: normalized jamming density of the su- 
perdisks, p.j(p)/pj{ 1/2), as a function of the deformation pa- 
rameter p. Upper panel: the p-derivative of the normalized 
jamming density near the special point p = 1/2. The symbols 
are the results of our MC simulations, whereas the solid lines 
show the approximation ([2}. 



3 



For supcrdisks, s(p) = S(l/2)/S(p) is given by follow- 
ing relations, 
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(4) 

Near p = 1/2, the approximation ([2]) is a continuous 
function of p, but has a jump in the p-derivative. In 
fact, the jump in the derivative of pj(p)/pj(l/2), given 
by our exclusion area approximation, is in a reasonable 
agreement with the result obtained by MC simulations 
presented in Fig. [3] The numerical values of the right 
and left p-derivatives of pj(p)/ pj(l/2) at p = 1/2 can be 



approximated by 

,. ds + (p) 
lim : 

p=i/2 dp 



lim 

0=1/2 



ds-(p) 
dp 



(5) 



In summary, in this work we demonstrated by numeri- 
cal MC simulations, as well as by approximate exclusion- 
area arguments, that the point at which the shape of the 
superdisks changes from concave to convex, is special not 
only in the geometric closed-packing properties, but also 
in the jammed-statc properties in RSA. Future work will 
be focused on the dynamical simulations, to explore the 
approach to the jammed state, as well as on studies of 
unoricntcd supcrdisks. 
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correspondence, and acknowledge support of this re- 
search by the NSF under grant DMR-0509104. 
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